Dynamical Dzyaloshinsky-Moriya Interaction in KCuFa 
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The spin dynamics of the prototypical quasi one-dimensional antiferromagnetic Heisenberg spin 
S = 1/2 chain KCuFs is investigated by electron spin resonance spectroscopy. Our analysis shows 
that the peculiarities of the spin dynamics require a new dynamical form of the antisymmetric 
anisotropic spin-spin interaction. This dynamical Dzyaloshinsky-Moriya interaction is related to 
strong oscillations of the bridging fluorine ions perpendicular to the crystallographic c axis. This 
new mechanism allows to resolve consistently the controversies in observation of the magnetic and 
structural properties of this orbitally ordered perovskite compound. 
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The pervoskite system KCuFa is presumably one of 
the best reahzations of an one-dimensional (ID) anti- 
ferromagnetic (AFM) Heisenberg chain as a direct con- 
sequence of the orbital ordering in this compound [l|. 
Recently, detailed inelastic neutron scattering investiga- 
tions have shown that fingerprints of a ID Luttinger liq- 
uid, namely the spinon excitation continuum, exist up to 
temperatures of about 200 K j^. Moreover, these specif- 
ically ID features are still observable below the onset of 
antiferromagnetic ordering at Tn = 39 K, indicating the 
strong quantum nature of magnetism in KCuFs. 

In spite of its paradigmatic status, the driving forces 
of the orbital and magnetic structure arc far from being 
understood at the moment. Even the crystal structure 
originally assigned to be tetragonal D^ (see Fig. [J) was 
claimed [Sj] to be orthorhombic D|. Electron spin res- 
onance (ESR) measurements [j, |5| played a key role in 
triggering these investigations: they suggested the ex- 
istence of the antisymmetric anisotropic exchange, usu- 
ally termed as Dzyaloshinsky-Moriya (DM) interaction, 
thereby questioning the crystal symmetry determined 
earlier. But the orthorhombic distortion deduced from 
x-ray diffraction [3| is not consistent with recent NQR 
data [6| and other experimental and theoretical findings 

In this Letter we reanalyze all static sources of the ESR 
line broadening on the microscopical level and show that 
the discrepancies in the understanding of this paradigm 
compound can be resolved by introducing a dynamical 
Dzyaloshinsky-Moriya (dDM) interaction. 

Single crystals of KCuFa have been grown as described 
in Ref. Isl. The ESR experiments were carried out with 
a Bruker ELEXSYS E500 CW-spectrometer at frequen- 
cies of 9 and 34 GHz. For measurements at 90 and 
150 GHz a quasi-optical technique was used [9| . The ESR 



spectrum consists of a single exchange narrowed Lorentz 
line at resonance fields corresponding to gc ~ 2.15 and 
ga = 2.27 in accordance wi th p revious experiments at 

Figure [2] shows the 
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24 GHz [10| and 34 GHz 
temperature dependence of the ESR linewidth for dif- 
ferent frequencies. Above Tn the linewidth increases 
monotonously from 0.3 kOe to more than 3 kOe above 
room temperature with a pronounced anisotropy with 
respect to the c-axis (see Fig. [ll^b)), but without any 
anisotropy within the ab plane in agreement with the 
all over tetragonal crystal structure. For different fre- 
quencies the linewidth data follow approximately the 
same temperature dependence with comparable absolute 
values. It can be described phenomenologically by an 
Arrhenius law Aif ex exp(— A/T) with an energy gap 
A = 114K (see Fig. O. 

Using conventional estimates 
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FIG. 1: (color online) (a): Schematic orbital ordering in the 
ground state of type-a KCuFa. The big (small) spheres denote 
Cu (F) ions. The K ions, placed between the CuFe octahedra, 
are not shown, (b): Two most relevant exchange paths of the 
ring-like SAE between the excited states of Cu ions. 
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FIG. 2: (color online) (a) Temperature dependence of the 
linewidth at different frequencies. The blue dashed line 
(34 GHz, H\\c) represents a fit to an Arrhenius law. (b) An- 
gular dependence of AH at T — 300 K. The fit is described 
in the text. 



sources of the ESR Unewidth, it was concluded that the 
relaxation contribution of the antisymmetric DM inter- 
action is several orders of magnitude larger than all other 
mechanisms in KCuFa [J|. However, the above men- 
tioned inconsistencies concerning the direction of the DM 
vector [3|, Ij, |5| and a general theoretical analysis of the 

1/2 chains 



spin relaxation in S ~ 1/2 chains [131 . |14| | questioned 
the validity of the original approach. Moreover, the mi- 
croscopic analysis of the exchange paths in several one- 
dimensional magnets [l5|, |l6|, [iTI, ll8| showed that the 
symmetric anisotropic exchange (SAE) may strongly ex- 
ceed the conventional estimate and thus dominate the 
spin relaxation. Therefore, it is necessary to reanalyze 
microscopically the magnitude of SAE and DM interac- 
tions in KCuFa. 

Let us reconsider the structure and the orbital ordering 
in KCuFs. The crystal has a pseudo-perovskite structure 
with a = 5.855 A, c = 7.852 A. Each Cu2+ ion (elec- 
tronic configuration 3(fi, spin 5 = 1/2) is surrounded by 
a tetragonally deformed fiuorine octahedron. The princi- 
pal axes of the octahedra are oriented alternating along 
the a and b axes, leading to the orbital ordering of the 



ground \y — z ) and 



) states of a hole on Cu ions 



for the anisotropic exchange and will be neglected in the 
following analysis. 

The crystal structure and the orbital ordering in 
KCuFa allow for SAE interaction W^^^ = 5^ ■ d''^^^ ■ S^ 
({a,/9} = {x,y,z}, a,b are the interacting ions) only 
along the crystallographic c = z axis via the ring-like 
processes described in Ref. |l7|. From the microscopical 



point of view, the exchange process consists of the ex- 
citation of a spin via spin-orbit (SO) coupling into an 
excited state of the starting ion a, the following hopping 
into an excited orbital state on the neighboring ion b 
and its transfer back into the initial state via the ground 
orbital state of the ion b. Fig. [T] shows that the over- 
lap between the ground and the relevant excited orbital 
states (|a;z), \yz)) via F~ p orbitals is of a and tt type, 
respectively. Denoting the respective hopping integrals 
by t'^ and t,r, we can estimate the exchange constant of 
SAE like in Ref. |I3 as 
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(see Fig. [Ua)). The existing small contribution of the 



where A is the SO coupling constant, Ac/ and Aab de- 
note the crystal-field splitting and the charge-transfer en- 
ergy, respectively. Here we used X„,/Ar,f ~ 0.05 [15[ and 
AUt'JAab « 4t;VAa6 « 380 K [3. 

The obtained value of SAE turns out to be strongly en- 
hanced as compared to the conventional estimate [a, ll9| , 
but the resulting ESR line broadening with AH ~ 10^ Oe 
cannot explain the huge ESR linewidth AH ^ 3 kOe ob- 
served in KCUF3. Note that in all compounds where 
SAE plays a dominant role [l5|, [la, lUl the character- 
istic linewidth at T ss J/k^ is about hundreds of Oe, 
i. e. an order of magnitude smaller than in the case of 
KCuFa. This large discrepancy indicates the existence 
of an additional source of line broadening in KCUF3. In- 
deed, Oshikawa and Affleck [i3| supposed another type 
of spin relaxation to be dominant in this system, because 
KCuFa does not obey the universal scaling behavior of 
AH{T/J) typical for spin relaxation via SAE. Therefore, 
we will turn now to the analysis of the second possible 
source of the ESR line broadening in KCuFa - the DM 
interaction: 

The existence of a static DM interaction is not allowed 
in the originally proposed tetragonal crystal structure [^ . 
However, AFM resonance at T = 4.2 K yields nonvanish- 
ing DM vectors d || [100] and [010] [1]. The reexamina- 
tion of the room temperature crystal structure by x-ray 
diffraction [3|] showed a large standard deviation of the 
fluorine ions placed on the magnetic chains away from 
the c axis. This was interpreted as a sign of a consid- 
erable static displacement of these ions from the c axis, 
leading to an orthorhombic distortion allowing for a non- 
vanishing DM vector directed along the [110] and equiva- 
lent directions. This kind of displacement in a perovskite 
structure usually means a rotation of a CuFg octahedron 



as a whole [20J , but no displacement of the F~ ions from 
the {xy) plane has been detected in Ref. |3|. This fact 
was emphasized by Binggeli et al. [7| who found in their 
LDA+U calculations that this structure relaxes to the 
conventional tetragonal structure of KCuFs. This indi- 
cates that the displacements reported in Ref. y are rather 
dynamic, initiated by thermal fluctuations. This assump- 
tion is corroborated by the measurements of the thermal 
displacement coefficients of fluorine ions [2]| . They were 
found to be very high and exceed strongly those of e. g. 
oxygen atoms in an oxide with the similar crystal struc- 
ture LaMnOa 2^- Such thermal motions may appear 



as static on the frequency scale of x-ray experiments 
/ ~ 10^'^ Hz. In the following we develop a model of 
dynamical anisotropic exchange, which accounts for the 
thermal vibration of an intermediate diamagnetic ion and 
allows to explain both the magnitude and the anisotropy 
of the ESR Hue in KCuFg. 

We start from the general expression for the antisym- 
metric exchange interaction TIdm = d^"'') • [Sa x Sf,]. In 
the fifth order of perturbation we have derived 
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where j = {x, y, z}, tap are the effective hopping integrals 
between the states |a), |/3) and a sum over all states of 
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FIG. 3: (color online) Possible paths for DM interaction be- 
tween two sites a (with ground state 77, exited state (") and 
b (^ and (p, respectively). Solid arrows correspond to the ef- 
fective hopping integrals, dashed arrows indicate the matrix 
elements of the spin-orbit coupling. Numeration indicates one 
of the possible order of the matrix elements in Eq. ((2]). Frame 
(I) displays additionally orbital configurations corresponding 
to this exchange process. 



an intermediate ion (|k), \p)) is implied. These virtual 
hopping processes are displayed schematically in Fig. [3l 
where e. g. frame (I) corresponds to the first term of 
Eq. (21). Note, that expression ^ differs strongly from 
the one suggested in Ref. |23|, where the position of the 
bridging ion, crucial in our case, was not discussed. 

It is well known [4|, |23| , that the static DM interaction 
is not allowed in the crystal structure of KCuFs with 
the space group Z)]^, but now we go beyond the static 
configuration. As illustrated in Fig. [31 the intermedi- 
ate ion c may be temporarily shifted from its equilib- 
rium position due to rotation or tilting phonon modes. 
If the average time of this displacement tf ~ 2Tr /lof is 
large compared to the characteristic time of electron ex- 
change Tab ~ Ti/tab, several exchange processes shown in 
Fig. [21 can occur via the displaced ion giving rise to an- 
tisymmetric exchange. Usually tab ~ 0.5 eV what allows 
to estimate Tab ~ lO^^rp using the phonon frequency 
ujp/2tt '^ 70 cm~^ as determined from Raman scattering 
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experiments 

The relevant pattern of the F^ displacements away 
from the c axis is taken in accord to x-ray diffraction 
experiments [3[. It can be related to the orbital or- 
der in KCuFa as follows: In Fig. [3l we show the geom- 
etry of exchange along with the ground state orbitals 
Iv) = \y^ - z'^)a and 1^) = |x^ - z^);,. There are four 
short and two long Cu-F bonds in a CuFg octahedron. 
A displacement of an intermediate F~ ion along the yc 
axis (as shown in the Fig. [3]) corresponds to the rotation 
of the Cu \x'^ — z'^)b orbital and fluorine ion around the 
Xb axis and does not destroy short bonds. Note that in 
this case the |y^ — z^)a orbital is fixed by three station- 
ary short bonds and is not rotated in accordance with 
the results of the structural analysis [3| which does not 
reveal any deviations of F~ ions from the (ab) plane. 
The second equivalent possibility, a displacement of the 
F~ ion along the Xc axis with the adiabatic rotation of 
the |y^ — z^)q around the ya axis, is not displayed in 
the Fig. [3l but will be taken into account. Substituting 
the excited state \ip) = \yz)a {{dyz\lx\dy'^-z^)a = 2i), the 
Eq. ([2]) gives only one nonzero component of d 
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Xa {2t',+t^){t',+t^)t'^ 



[-t'^cos^{29)+t^sm'^{2e)]cos{29)sm{2e). 



(3) 



Note that the factor sm{29) can be rewritten as a modu- 
lus of the vector product [uac x Ubd , where the unit vec- 
tors Uac and iibc connect the spins a and b with the bridg- 
ing ion c, respectively. Therefore, the rule d || [uac x ribc] 
2t 



I for the static case is preserved 



suggested in Refs. |25 
in our case, too. 

The value of d is related to the displacement of a F^ 
ion from equilibrium position sin(20)-i?cu-F, where i?cu-F 
is the distance between the Cu and F ions and can be ex- 
pressed via phonon operators. Thus we arrive to a new 



dynamical form of DM interaction (dDM) containing si- 
multaneously a phonon and two spin operators. A closely 
related form of spin-lattice interaction was introduced by 
Kochelaev [27| (see Eq. (5) in Ref. |28[ ) in the context of 
explanation of the isotope effect on the ESR linewidth in 
lightly doped La2Cu04. In our case, however, we do not 
need to recall the conventional spin-phonon interaction 
(coupling parameter G in Refs. 1271 . 1281 ) at all, which was 
a crucial ingredient in that approach. 

In the following, we will use this dynamical Hamil- 
tonian to calculate the ESR linewidth by the method 
of moments Ai7 oc M2{(P)/J [H, [H, [l3|, where upon 
the thermal averaging over the phonon variables will be 
performed. Applying the Einstein model of vibration 
we get (sin^(26')i?^^_p) = ;^coth(?iwF/2A:BT), where 
m denotes the mass of an oscillating fluorine ion and 
wf is the frequency of the rotating mode. The quantity 
(sin(26')i?cu-F) is naturally equal to zero and the static 
DM interaction does not exist. To compare with the re- 
sults of x-ray analysis, it is useful to introduce the root- 
mean-square deviation of the F~ ion away from the c 

axis Ai?^ = i/(sin^(26')i?^j^_p). Using the typical ro- 
tation frequency of CuFg octaliedra wp ^ 70 cm~^ [2J| 
one can estimate Ai?j_ w 0.07 A at T = 100 K, what 
coincides with Ref. y. Therefore, we conclude that these 
displacements are rather dynamic than static at T > T^. 
This conclusion is supported by the latest NQR study [6[ 
which docs not give any evidence for a static dislocation 
of F~ ions from the z axis at 77 K as well. To estimate 
dx we will use the values given in Eq. ([3]). Recalling 
the second possible direction of motion of intermediate 
F^ ions, we derive finally 
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5.6 K. 



(4) 



Note that this estimate coincides very well with the value 
'^AFM ~ ^-13 K employed by Yamada and Kato [5| to 
describe the anisotropy of the AFM resonance line below 
Tn. Moreover, it yields the correct linewidth of about 
10'^ Oe in the high-temperature limit. 

The anisotropic interchain exchange can be neglected, 
because of the small constant of isotropic exchange Jj^ ^ 
J and the orthogonality of orbital states. The dynamical 
form of SAE would arise in the higher order of perturba- 
tion theory as compared to the dDM interaction and is 
expected to be negligibly small. The anisotropic Zecman 
effect does not contribute, because we do not observe 
any significant field dependence of the linewidth up to 
the frequency of 150 GHz (Fig. [2|). Therefore, the spin 
relaxation mechanisms described above, the static SAE 
(Eq. [T|) and the dDM interaction (Eq. [4]) between the Cu 
ions along the c axis are the only relevant sources of the 
ESR line broadening in KCuFa. Indeed, these two broad- 
ening mechanisms can explain both the magnitude and 
the anisotropy of the ESR linewidth. Figure [2I^b) shows 
a fit of the angular dependence of A_ff at room tempera- 



ture and emphasizes the contributions of SAE and dDM 
interaction to the linewidth, which reflects clearly the 
dominant role of the dDM interaction. This indicates 
that KCuFs, despite its paradigmatic status for orbital 
order, is governed by strong fluctuations in the lattice 
and orbital sector. 

In summary, we have introduced a dynamical 
Dzyaloshinsky-Moriya (dDM) interaction which allowed 
to explain the long standing mystery about the origin of 
the huge ESR linewidth in KCuFs at T > Tn in consis- 
tence with other experimental findings about the proper- 
ties of this quasi-lD system. Such a dDM interaction be- 
comes effective, if the characteristic time of the dynamic 
distortion resulting in a nonzero DM vector is large com- 
pared to the time scale of the exchange interaction and if 
the amplitude of these distortions is high. This is the case 
for low-lying optical modes with the tendency to soften 
to low temperatures. Therefore, such kind of interaction 
may be on general useful for understanding the spin dy- 
namics of other materials with soft-mode vibrations. 
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